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ABSTRACT 

Modified Newtonian dynamics (MOND) is an alternative theory of gravity that aims 
to explain large-scale dynamics without recourse to any form of dark matter. However 
the theory is incomplete, lacking a relativistic counterpart, and so makes no definite 
predictions about gravitational lensing. The most obvious form that MONDian lensing 
might take is that photons experience twice the deflection of massive particles moving 
at the speed of light, as in general relativity (GR). In such a theory there is no general 
thin-lens approximation (although one can be made for spherically-symmetric deflec- 
tors), but the three-dimensional acceleration of photons is in the same direction as the 
relativistic acceleration would be. In regimes where the deflector can reasonably be 
approximated as a single point-mass (specifically low-optical depth microlensing and 
weak galaxy-galaxy lensing), this naive formulation is consistent with observations. 
Forthcoming galaxy-galaxy lensing data and the possibility of cosmological microlens- 
ing have the potential to distinguish unambiguously between GR and MOND. Some 
tests can also be performed with extended deflectors, for example by using surface 
brightness measurements of lens galaxies to model quasar lenses, although the break- 
down of the thin-lens approximation allows an extra degree of freedom. Nonetheless, it 
seems unlikely that simple ellipsoidal galaxies can explain both constraints. Further, 
the low-density universe implied by MOND must be completely dominated by the 
cosmological constant (to fit microwave background observations), and such models 
are at odds with the low frequency of quasar lenses. These conflicts might be resolved 
by a fully consistent relativistic extension to MOND; the alternative is that MOND 
is not an accurate description of the universe. 

Key words: gravitational lensing - relativity - gravitation - dark matter - acceler- 
ation of particles. 



1 INTRODUCTION 

Modified Newtonian dynamics (MOND; Milgrom 1983a) is 
an alternative theory of gravity in which either inertia or 
the effective gravitational force is changed in the limit of 
low accelerations (as opposed to large distances). Taken in 
its purest form, a MONDian universe would contain no dark 
matter, and be dominated by the low density of baryons 
generated in primordial nucleosynthesis (Sanders 1998; Mc- 
Gaugh 1999). MOND makes definite and generally success- 
ful predictions about the dynamics and properties of galax- 
ies (Milgrom 1983b; Mateo 1998; McGaugh & de Blok 1998; 
Sanders 2000), groups (Milgrom 1998), clusters (Milgrom 
1983c; Sanders 1994) and large-scale structure (Milgrom 
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1997). However MOND is not a complete physical theory 
- in particular it has no relativistic counterpart (Milgrom 
1983a; Bekenstein & Milgrom 1984; Sanders 1997). Thus 
there are no concrete cosmological predictions, and it is not 
certain how photons couple to gravitational fields. 

Despite the lack of any relativistic theory underpinning 
MOND, it is possible to constrain its properties by invok- 
ing the appropriate Newtonian limits and using basic sym- 
metries. Felten (1984) and Sanders (1998) argued that the 
universe as a whole is in the Newtonian regime and so obeys 
the standard Friedmann equations^] and similarly that the 

t Thus the cosmology can be specified by the usual parameters: 
the normalised present day matter density, f2 mo , the similarly 
normalised cosmological constant, Ha , and the Hubble constant, 
H = 100ft km s" 1 Mpc -1 . 
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early universe is governed by general relativity (GR). Tak- 
ing Ho = 70 km s _1 Mpc -1 , and assuming a totally bary- 
onic universe, the nucleosynthesis constraints of Tytler et 
al. (2000) imply that fi mQ ~ 0.01. McGaugh (2000) showed 
that the CMB anisotropies measured by the MAXIMA-1 
(Hannay et al. 2001) and BOOMERanG (de Bernardis et 
al. 2000) experiments are consistent with such a low-density 
cosmological model provided that the universe is spatially 
flat, i.e., S7a = 1 — f2 mo — 0.99. Despite these successes, 
there are a number of qualitative difficulties with a MON- 
Dian cosmology (e.g. Scott et al. 2001), although some are 
merely indicative of the great conceptual differences between 
MOND and more conventional physics. 

There has been less focus on gravitational lensing within 
a MONDian framework, primarily as the deflection of light 
cannot be explained without recourse to a relativistic the- 
ory. Beckenstein & Sanders (1994) and Sanders (1997), in 
attempting to derive a relativistic extension of MOND, in- 
cluded gravitational lensing in their considerations, but a 
more empirical approach is adopted here. Mortlock & Turner 
(2001) used the galaxy-galaxy signal measured by Fischer et 
al. (2000) to constrain the deflection law of galaxies to be 
A(R) oc i? 01±0 1 over the range 10 kpc <> R <, 1 Mpc. This 
is consistent with the deflection law obtained by assuming 
that, as in GR, photons experience twice the deflection of 
massive particles moving at the speed of light (Qin, Wu & 
Zhu 1995). Thus the asymptotic lensing effect of a MONDian 
point-mass matches that of an isothermal sphere, which in 
turn is consistent with most observations of lensing by galax- 
ies (e.g. Brainerd, Blandford & Smail 1996; Kochanek 1996). 

Encouraged by the qualitative agreement between data 
and observations, the lensing formalism of Qin et al. (1995) 
is explored in more detail here. In Section ^| a number 
of general results are derived and then applied to simple 
lens models. These are then compared to observations of 
galaxy-galaxy lensing (Section^), microlensing (Section^), 
multiply-imaged quasars (Section ^ and lensing by clusters 
(Section ^|). The conclusions are summarised in Section ^ 



where ao is the critical acceleration, above which the dy- 
namics is close to Newtonian. The exact form of fyi(x) is 
not known, but /m(x) = 1 if x 3> 1 and /m(x) = x if x <C 1. 
No major MONDian results are sensitive to the form of this 
function in the regime x ~ 1, and so 



/m(z) = 



x, if x < 1, 



1, if x > 1 



(3) 



is used here for mathematical simplicity. It is convenient to 
invert equation (^) so that the acceleration can be expressed 
as a function of known quantities. Defining a new function 
f-U(x) with the same asymptotic limits as /m(x), MOND 
may also be defined by 



a(r) = fZ, 



-1/2 



a N (r) 



Oo 



a N (r) 



(4) 



For the particular choice of /m(x) given in equation ([J), 
/m(z) = Im{x). 

There are some difficulties in understanding equa- 
tions (|^) and (^) for complex systems (e.g. Milgrom 1983a; 
Scott et al. 2001), but their qualitative intrepretation for 
a single particle is reasonably straightforward. For accel- 
erations much greater than ao Newtonian dynamics re- 
mains valid, as a(r) = a^(r). But if an(r) <C ao then 
a(r) = [aoa^(r)] lll2 a-^(r) and a test particle would feel an 
acceleration in the expected direction, but of a greater mag- 
nitude. The value of the critical acceleration has been deter- 
mined from galaxy rotation curves as ao = 1.2 ±0.1 x 10~ 10 
m s -2 (Milgrom 1983b; Begeman, Broeils & Sanders 1991; 
McGaugh & de Blok 1998), and, modulo the uncertainty in 
the form of 'fM(x), MOND has no other free parameters. 
Equations ^ and (^) are sufficient to calculate the trajec- 
tory of a massive particle (Section 2.1), from which a natural 



formalism for the def lection of photons can be extrapolated 
(Sections O andO). 



2 MOND 

MOND is most generally a modification of inertia (Milgrom 
1983a), but, for gravitationally-dominated systems, can be 
more intuitively treated as an increase in the gravitational 
acceleration felt by test particles. The Newtonian accelera- 
tion due to a (static) mass distribution p(r) is 



a N (r) 



-G 



p(r) 



(1) 



where G is Newton's gravitational constant and r is posi- 
tionjfj The MONDian acceleration is related to this by (Mil- 
grom 1983a) 



a(r)fu 



a(r) 



ao 



(2) 



t Three-dimensional vectors are denoted by lower case (e.g. r); 
two-dimensional vectors - projected onto the sky - are denoted 
by upper case (e.g. Ft). 



2.1 Deflection of massive particles 

Within the framework of MONDian (or Newtonian) dynam- 
ics it is possible to calculate the trajectory of an arbitrary 
particle. If, however, the particle is sufficiently light that 
it does not disturb the external mass distribution, and its 
speed is such that its path is nearly linear, its deflection can 
be approximated by integrating the acceleration along the 
unperturbed trajectory (i.e. a straight line). Without loss of 
generality, this path can be defined as being parallel to the 
z-axis of a Euclidean coordinate system, with R = {x, y) the 
two-dimensional impact parameter (relative to some refer- 
ence point) in the plane of the sky. Thus the deflection angle, 
A — (A x , A y ), of a particle moving with speed v in the z- 
direction is given by 



Ai(R) 



a z [r(t)]dt 

-i/2 / a n[r(t)] 
ao 



a m [r(t)] dt, 



(5) 
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where i € {x,y}, and a(r) is denned in equations ([l]) and 
(^). For a minimally-deflected particle, r(t) = [x, y, v(t—to)], 
where to is an arbitrary reference time. 

2.2 Deflection of photons 

In the standard theory of gravitational lensing (e.g. Sch- 
neider, Ehlers & Falco 1992), the deflection of photons in 
weak, static gravitational fields can be calculated from GR, 
and is simply twice the deflection experienced by a mas- 
sive particle moving at the speed of light, c. This has been 
confirmed observationally for light grazing the Sun (Dyson, 
Eddington & Davidson 1920; Robertson & Carter 1984) and 
this simple relationship must also be true of MONDian lens- 
ing in the Newtonian limit. There is no strong evidence that 
this relationship holds in the deep MONDian limit, but such 
a model is consistent with galaxy-galaxy lensing (Mortlock 
& Turner 2001), and this hypothesis is adopted henceforth. 
The deflection angle of photons can thus be read from equa- 
tion (^) as 

MR) = § f /m" 1/2 {^i}a Ni [r(t)]d*, (6) 

where again i £ {as, y} and, for a photon, r(t) = [x,y,c(t — 
to)]- 

2.2.1 The thin-lens approximation 

In the limit ao — > 0, the standard relativistic deflection for- 
mula can be reproduced by using equation (|l|) to obtain 



A Ni (R) = - 

c 



ctNi[r(t)] dt. 



(7) 



In the limit of small deflection angles (relative to the line- 
of-sight extent of the lens) the thin-lens approximation (e.g. 
Schneider et al. 1992) can be invoked. Equation (R) can thus 
be reduced to 



A N (R) 



4GS(fl') R' R 
c 2 \R' R\ 2 



d 3 R', (8) 



where the surface density of the lens is 
= / p(x,y,z)dz. 



(9) 



The application of the thin-lens approximation allows man 
lensing problems to be greatly simplied, and equation 
is one of the main results of standard gravitational lensing 
theory. 

As given in equation (^) the MONDian deflection angle 
cannot be simplified in general. In particular, there is no use- 
ful thin-lens approximation: equation (^j) cannot be simpli- 
fied in the same way that equation ([?]) can. More qualitative 
arguments imply that the same is lik ely to be true of any 
MONDian lensing theory (see Section 12.3.2). Further, con- 



stant mass-to-light ratio models of real gravitational lenses 
cannot reproduce the observed image configurations (e.g. 
AbdelSalam, Saha & Williams 1998; Section [Tl]) ; this would 
almost certainly invalidate MOND if the thin-lens approxi- 
mation were valid. The ambiguity in the line-of-sight mass 
distribution of an observed deflector is inconvenient (as is 
the extra integral required in the calculation of deflection 



angles), but could be thought of as analogous to the a pri- 
ori unconstrained dark matter distributions that dominate 
conventional gravitational lensing studies. 

Despite the lack of a general thin-lens approximation, 
some simplifications of equation (^) are possible. In many 
cases the deflection angle can be decomposed into two MON- 
Dian segments (when the photon is far away from the lens) 
and a single Newtonian segment (the major deflection near 
the lens). If the deflector is sufficiently localised, the latter 
regime can be handled in the usual way, using the thin-lens 
formalism, and the MONDian contribution away from the 
lens ma y be ap proximated by the simple point-mass formula 
(Section 2.3.1), the deflector being 'unresolved' at such great 
distances. 



2.2.2 The lens equation 

Gravitational lensing can be thought of as a mapping be- 
tween the source plane (an imaginary surface, perpendicular 
to the line-of-sight, on which a source is located) and the sky 
or image plane. The lens equation relates the angular posi- 
tion of a source, (3, to the angular position(s) of its image(s), 
0, and can be written as (cf. Schneider et al. 1992) 



j3 = + a(0), 



where 



fids 



a(0) = -^A(d od G) 



(10) 



(11) 



In these definitions, d d, d os and dds are the angular 
diameter distances between observer and deflector, observer 
and source, and deflector and source, respectively. It is not 
entirely clear how they vary with redshift in a MONDian 
cosmology. Even if the universe obeys the Friedmann equa- 
tion, as implied by the results of Sanders (1998), the sepa- 
ration of nearby photons can only be determined by a rela- 
tivistic extension of MOND, although the standard distance 
measures should be reproduced in the limit of large angu- 
lar separation (cf. Linder 1998). Where values of d a d, d os 
and dds are required, they are calculated for the appropriate 
low-density Friedmann model using standard formulae (e.g. 
Carroll, Press & Turner 1992). This is probably consistent 
with the lensing theory described above, but is no more than 
a reasonable assumption. 

The lens equation allows the calculation of image posi- 
tions, magnifications and distortions. The magnification of 
a single image is (Schneider et al. 1992) 



ti(0) = 



d 2 /3 



d0 2 



(12) 



and the total magnification of a source is then 

Mtot(j9) =X>[0<G9)], (13) 

i=l 

where the sum is over the Ni images formed, and 0i(j3) is 
the position of the ith image. 



2.2.3 Spherical symmetry 

If the deflector is spherically symmetric, then p(r) = p(r) 
and equation (H) can be written as 
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MR) = 



/m 



where 



2G 



R 



(14) 



-1/2 



GM(< VR 2 + z 2 
(R 2 + z 2 ) a 



M (< VR 2 + zA 



(R 2 + z 



,2^/2 



dz, 



M(<r) = / 4nr' p(r')dr'. (15) 
Jo 

It is possible to write down an analogue ol the usual lens 
equation by expressing M(< r) in terms of the projected 
surface density. Using equation ( |l5|) and expressing p(r) in 
terms of an Abel integral (e.g. Binney & Tremaine 1987) , 



M{< r) 



4irr 



1 dE 



7T dR 



VR? 



.dRdr\ 



R? 



R^dR 



dR 



(16) 



where the order of integration is reversed, making the inner 
integral trivial (cf. Kovner 1987). This formula is not par- 
ticularly useful for unbou nded mass distributions like the 
isothermal sphere (Section 2.3.3), as its convergence proper- 
ties are poor. However in the case of MOND galactic mass 
distributions fall off as fast as their luminosity density - 
exponentially, in general - and so equation ( |l6| ) is well be- 
haved, as well as relating the observed surface brightness of 
a deflector directly to its lensing properties. 

The combination of equations (jlj) and ( |lrj ) are sugges- 
tive of an alternative formalism for spherical lenses. MON- 
Dian lensing can be cast in terms of the conventional formal- 
ism of Schneider et al. (1992) if an 'effective' surface density 
is defined. Specifying it by either Pn(t") or En(-R) (as it is 
only relevant in the symmetric case) , it is the density profile 
that would, under the assumption of Newtonian mechanics, 
produce both the same rotation curve and the same deflec- 
tion law as the model profile produces in MOND. Given that 
a point-mass has an asymptotically constant rotation speed 
of v c = (GMa ) 1/4 in MOND (Milgrom 1983b), the addition 
of an isothermal component is suggested, and defining 



pn{r) = p{r) + ■ 



(17) 



47rG l/2 r 2 

is a reasonable approximation to many MONDian deflec- 
tors of true mass M. The same rotation curves and deflec- 
tion laws would result from defining fyi(%) = 1 + [1 + (1 + 
4x) 1/2 ]/(2x). 

In the limit of /m(x) — > 1 and ao — > 0, equation ( |l4"| ) 
reduces to the conventional result that (e.g. Schneider et al. 
1992) 



A N (R) 
where 



4G Af (< R) 
' c 2 R 2 



R, 



M{< R) = 



R)= [ 2tt_R'E(J?') dR'. 
Jo 



(18) 



(19) 



s 




1000 



Figure 1. The deflection law of a point-mass in MOND (solid 
line) and GR (dashed line). The dotted line shows the approxi- 
mation to the MONDian result obtained by calculating the New- 
tonian deflection caused by an iso thermal sphere superimposed 
on a point-mass (see Section 2.2.3). 



2.3 Lens models 

Having developed a plausible formalism for gravitational 
lensing within the framework of MOND, more specific re- 
sults can be derived for some simple mass distributions. 



2. 3. 1 Point-mass 

For a point lens of mass M, p(R) = MS 3 (R) and M(< R) = 
M. Using the piecewise definition of fyi(x), equation ( |14[ ) 
splits up into several trivial integrals, and the deflection law 
for a point-mass is given by 



A(R) = 



(20) 



4GM 
'~c r R 



B 2 



1 arctan 



- 1 



if R < r M , 



—Am if R > t"m, 

where tm = (GM/ao) 1 ^ 2 is the distance from a point-mass 
at which the physics changes from the Newtonian to the 
MONDian regime and A M = 2-K(GMa ) 1/2 /c 2 is the asymp- 
totic MONDian deflection angle. This is compared to the 
standard point-mass deflection angle (as inferred from the 
Schwarzschild metric) in Fig. |l|. For large impact parame- 
ters (R <; f m) the photon experiences an acceleration of less 
than ao, and so the deflection angle is independent of impact 
parameter, in much the same way that the rotation speed 
would be. For R ^ tm, the light path is in the MONDian 
regime for the most part (the second term in the above for- 
mula), but the deflection is dominated by the Newtonian 
force close to the lens (the first term, which approaches the 
Schwarzschild bending angle for small R). 

The next step in analysing the point-mass lens is to 
convert equation ( po| ) into angular units, in order to write 
down the lens equation. The Schwarschild lens is completely 
characterised by its Einstein angle, 
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4GM d ds 



d od d a 



(21) 



(22) 



but two other angular scales are uselul here: 
8m = — = / GM 

beyond which the deflection angle becomes constant; and 
2ny/GMa d ds tt 6>| 



do 



2 9k 



(23) 



the value of the deflection angle in this regime. Combining 
the above definitions with equation (jllj) , equation ( pp| ) re- 
duces to 



*(<?) = 



(24) 



S2 



• arctan 



— qm, if 6 > 6 M , 

with a(0wi) ~ -qm. From equation (|l^) the magnification 
of a single image is then 



(25) 



2 6* M + ^eV^m" ^ 2 



1 arctan 



){< 



L 61/(61 ■ 



CUM) 



-e 2 

< 9m, 
> Om, 



whereas for a Newtonian point-mass it is simply /xn(#) = 

\e A /{e A -ei)\. 

Inverting the lens equation [equation (11 01) ] and u sing 
the magnification formula given in equation(|l3|) and ( |25[ ) 
allows fi to t to be calculated. This is shown in Fig. ^, with 
the x-axis normalised to match the Newtonian Einstein ra- 
dius. In all cases /itot(/3) matches the relativistic form for 
(3 ^ 9m; beyond that MOND always results in a greater 
magnification. The shapes of these curves are determined 
by the mass-independent ratio of the Einstein angle to the 
characteristic MONDian angular scale. From equations (J2TJ) 
and (E2) this ratio is given by 



0e_ 

0m 



4a d od d d 



do 



<, 0.006 



1 Mpc 



1/2 



: 0.352, 



1/2 



(26) 



where the upper limit is calculated by assuming that dod — 
dds — dos/2, and the last expression is valid only in the 
local universe, where d os = cz s /Ho- This ratio is shown in 
Fig. |i| its increase with z B is quite apparent, and it is also 
important that it does not vary greatly with z d unless either 
ii d or dds is close to zero. Within the Galaxy MONDian 
lensing effects should be minimal, with 6e/8m Si 0.001, but 
they could become very important on cosmological scales, a 
point which is explored in more detail in Section ^. 

2.3.2 Thin rod 

In the standard theory of gravitational lensing, the deflection 
properties of most isolated mass distributions depend only 




1000 



1000 



Figure 2. The total magnification of a source as a function of its 
angular separation, (3, from the optical axis of a MONDian point- 
mass deflector. More extreme effects are emphasised in (a) and 
weaker effects shown in (b). The solid curves show increasing val- 
ues of 0e/$M; as labelled; the dashed curves show the Newtonian 
result (which could also be obtained by taking #e/#m = 0)- Note 
that 0e is the Einstein angle of the Newtonian point-mass; the 
x-axis could also be normalised by the actual Einstein angle of 
the lens, which increases as 9m decreases. This approach is taken 
in Fig. 



on their projected surface density. As argued in Section 2.2 
the thin- lens approximation is unlikely to hold in MOND, 
and this can be illustrated by considering the lensing effect 
of a thin rod oriented along the line-of-sight. Its density is 
given by 



0, 



if z < -1/2, 



p(r) = { M/IS 2 {R), if -l/2<z<l/2 



(27) 



0, 



if z > 1/2, 



where M is the rod's mass and I its length. Inserting this 
definition into equation ([]) and evaluating the required in- 
tegrals gives the Newtonian acceleration as 
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Figure 3. The ratio of the Einstein angle of a point- mass, 6^, to 
its characteristic MONDian angular scale, 9m i as a function of of 
source redshift, z B (a) and deflector redshift, z^ (b). In (a) the de- 
flector is assumed to be halfway between observer and source, and 
in (b) results are shown for several source redshifts, as indicated. 
Two cosmological models are shown: f2 mo = 0.01 and S7a = 0.99 
(solid lines); and Q mg = 0.01 and Ha = 0.0 (dashed lines); the 

1/2 

dotted line in (a) is given by #e/$m = 0.352 s , which is a good 
approximation in the local universe. Scales corresponding to the 
Galactic centre, the Magellanic Clouds, and the local group are 
marked, with the 'effective redshift' given by z B = Ho/cd OB . 



On, 



ON.= - 



GM x 

'~rn? 

GM y 
GM 



z + 1/2 



z - 1/2 



sjB? + (z + 1/2) 2 yjB? + (z - l/2f 
z + l/2 z-l/2 



I 



^R 2 + (z + l/2) 2 y/R 2 + (z- 1/2) 
1 1 



y/R 2 + (z- 1/2) 2 y/IP + (z + 1/2) 2 



(28) 



where, as before, R 2 — x' 2 + y 2 . The MONDian acceleration 
can be calculated directly from this using equation (^), but 
the results are cumbersome and so not presented in full. 




100 1000 



Figure 4. The equatorial acceleration produced by a thin rod 
(oriented along the z-axis) in MOND (solid lines) and GR (dashed 
lines). Here is the radius at which the acceleration becomes 
MONDian for the corresponding point-mass. Curves are shown 
for several lengths: 0; O.lrM; fM! lOrnj; and 100rM, with the 
acceleration decreasing (at least for r < rjyi) as length increases. 



Some idea of the lensing properties of such an object can 
be gained by looking at the acceleration in the x-y plane. 
Applying the definitions given in Section ^| yields an 'equa- 
torial' MONDian acceleration (directed towards the origin) 
of magnitude 



^£[i? 2 + (//2) 2 ] 



21-V2 



a{R)\ 



where 



^ \R 2 + (l/2) 2 ]- 1/4 , 



I 1 

2^2 



1 + 4 



T/2 



1/2 



if R < R M , 

if R > R M , 

(29) 



(30) 



is the equatorial radius at which the acceleration equals ao, 
and tm is the MONDian radius of the equivalent point-mass. 
A family of these curves is shown in Fig. with the New- 
tonian accelerations also included for comparison. Close to 
the rod a(R) oc R~ , the standard Newtonian result, and for 
very large impact parameters a(R) oc R~ , the rod acting 
like a MONDian point-mass. In most cases there is an inter- 
mediate region in which the acceleration is sufficiently small 
to feel a MONDian boost but the rod's finite length is also 
important. The combination of these two effects results in a 
a(R) oc i? _1//2 dependence, and any particle passing through 
this region would feel an unusually large gravitational pull. 

It is also revealing to consider how the acceleration 
varies along the z-axis. Taking x = y = in equation ( |28[ ) 
and again applying MONDian dynamics (outside the length 
of the rod) yields 



a(0,0, 



(31) 
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GM 
I 



-1/2 \z\+l/2 



if 1/2 < \z\ < ZM, 



\z\-l/2 



\z\+l/2 



1/2 



z, if 



\z\ > Z M , 



where 



zm 



1 + 



1/2 



(32) 



Note that if I ^> tm the acceleration becomes MONDian for 
\z\ 1/2, i.e. anywhere outside the rod. 

The integrals required to obtain the deflection angle 
must be calculated numerically, and the results are shown 
in Fig. If R <C Rm then the deflection is essentially 
Newtonian [i.e. A{R) — > AGM/(c 2 R)] as the deflection is 
dominated by the portion of the photon path alongside the 
rod. The opposite extreme is if R ^> /, in which case the 
rod acts like a point-mass, being 'unresolved' by the pho- 
ton. Thus A(R) -> 2n(GMa ) 1/2 /c 2 = A M in this regime, 

provided that the rod is of finite length^ A corollary of 
these two arguments is that the thin-lens approximation is 
valid in MOND if the line-of-sight extension of the deflector 
is less than tm, as, for all impact parameters, the deflec- 
tion is either Newtonian or effectively point-like. However if 
i iJ tm then A(R) oc RT 1 ^ 2 for Rm & R & I. In this case 
the photon's deflection is dominated by a segment of length 
~ I during which it feels a nearly perpendicular accelera- 
tion of a ~ [2GMa /(Rl)] 1/2 . Thus the deflection angle is 
A(R) ~ 2/c 2 (2GMa l/R) 1/2 in this region, and is propor- 
tional to the line-of-sight extent of the deflector. 

The differences between the curves in Fig. |H| give some 
idea of the the validity of the thin-lens approximation for 
an isolated deflector. For example, a typical elliptical galaxy 
with its major axis along the line-of-sight would have an 
effective / of ~ 10 kpc or more (Section ^) and tm — 2 
kpc, which would imply that the thin-lens approximation is 
probably reasonable in this situation. 



2.3.3 Isothermal sphere 

The isothermal sphere (e.g. Binney & Tremaine 1987) is 
a useful and simple galaxy model within the dark matter 
paradigm as it explains the flat rotation curves of spirals, the 
dynamics of ellipticals, and the lensing properties of both. 
Parameterised by its (unobservable) velocity dispersion, a, 
and its core radius, r c , its density is given by (e.g. Hinshaw 
& Krauss 1987) 



PN(r) 



(33) 



2ttG R 2 +r 2 ' 

Given the choice of /m(:c) described in Section |^, the 
gravitational properties of this model can be almost exactly 
replicated in MOND by the density profile 



p{r) = 



2ttG R?+r'i ' 



if r < r M , 



if r > r M , 



(34) 



3 An infinitely long rod has an infinite deflection angle, irrespec- 
tive of impact parameter. 



s 




1000 



Figure 5. The deflection law of a thin rod oriented along the 
line-of-sight in MOND (solid lines) and GR (dashed line). Here 
r'M is the radius at which the acceleration becomes MONDian 
for the corresponding point-mass, and Ayi is the asymptotic con- 
stant deflection angle. Curves are shown for several lengths: 0; rM 
(almost indistinguishable from 1 = 0); 10rM; 100rM; and IOOOm, 
with the deflection angle increasing with length. 



assuming r c < rij, which is certainly the case for real galax- 
ies. The image positions and magnifications produced by 
such a mass distribution could be obtained from the MON- 
Dian formalism (Section ^|), but there is no need to per- 
form these calculations explicitly, as its lensing properties 
match those of the conventional isothermal sphere, which 
has been studied in great detail (e.g. Hinshaw & Krauss 
1987; Kochanek 1996; Mortlock & Webster 2000). In partic- 
ular, standard calculations of lens statistics remain valid, a 



point explored in more detail in Section 5.2 



3 GALAXY-GALAXY LENSING 

Background galaxies are observed to be tangentially aligned 
around foreground galaxies due to the latter population's 
gravitational lensing effect. Termed galaxy-galaxy lensing, 
this technique has - under the assumption of GR - con- 
firmed the existence of extended isothermal haloes around 
all types of galaxies (e.g. Brainerd et al. 1996; Fischer et 
al. 2000). However it is potentially an even more power- 
ful probe of alternative gravity theories as the shear signal 
extends far enough beyond the visible extent of the deflec- 
tors that they can be treated as point-masses. The radial 
dependence of the shear signal is consistent with the MON- 
Dian lensing formalism described in Section ^| (Mortlock & 
Turner 2001), and the forthcoming Sloan Digital Sky Survey 
(SDSS; York et al. 2000) data should allow measurements 
out to physical separations of ~ 1 Mpc. Beyond this the in- 
fluence of secondary deflectors will start to dilute the signal, 
but, if MOND returns to a Newtonian regime at ultra- low 
accelerations (~ 10~ 12 m s -2 ), the resultant cut-off might 
be detectable. 

The full SDSS data-set will also facilitate a number 
of more general tests which should be able to discriminate 
between dark matter and alternative gravity theories un- 
ambiguously. These are discussed in detail by Mortlock & 
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Turner (2001), and so only summarised here. One implica- 
tion of MOND is that the shear signal around any subset of 
the foreground population should have the same functional 
form, within the errors; any departure from this would repre- 
sent strong evidence against such a model. A more powerful 
idea is to search for any deviation from azimuthal symmetry 
in the lensing signal (cf. Natarajan & Refregier 2000). A pos- 
itive detection would be difficult to reconcile with MOND 
(or any alternative gravity theory with basic symmetry prop- 
erties); a symmetric signal, on the other hand, would be 
in conflict dark matter-based halo formation models (e.g. 
Navarro, Frenk & White 1995, 1996). 



4 MICROLENSING 

Microlensing is the term used to describe the lensing action 
of (collections of) individual stars, and has been observed 
over a wide range of scales, from within the Galactic halo 
to redshifts of order unity. As the name suggests, the resul- 
tant image separations are much smaller than present day 
telescopes can resolve, and its only observable consequence 
is the change in magnification of a source caused by the 
relative motion of the deflector(s) across the line-of-sight. 

Given that the additive nature of MONDian lensing is 
unknown, this discussion of microlensing will be limited to 
isolated deflectors, primarily single stars. For a point-source 
of unlensed magnitude mo, the resultant light-curve takes 
the form (cf. Paczyhski 1986) 



m(t) 



mo 



(35) 



2.5 log < (Mm 




where /itot(/3) is given in Section 2.3.1, /3 m in//3o characterises 
the greatest alignment of deflector and source, to is the time 
required for the alignment to change by angle f3o , and i m i n is 
the time at which the greatest alignment occurs. The choice 
of /3o is somewhat arbitrary, but it is usual to define f3o — 9e , 
which implies that /itot(A)) = 1-34 in GR. This last relation 
is taken to define /3o here, which then implies that /3o > #e 
in general, with equality in the limiting case of 9m — > oo. 

Several MONDian microlensing light-curves are shown 
in Fig. []. The distinction between MOND and GR is clear if 
Qe/0m ii 1; otherwise the only difference between the light- 
curves is in the m(t) ~ mo wings. The ratio Be / 9 m de pends 
mainly on the observer- source distance (Section 2.3.1), and 
only approaches unity if z s <; 1 (Fig. ^). In physical terms, 
microlensing is useful as a probe of gravitational theories 
if the magnification is significant along geodesies passing 
beyond the Newtonian region of the deflector. This is no t th e 
case for microlensing within the local group (Section |4.l[ ), 
but is tr ue ( in a MONDian universe) on cosmological scales 
(Section^). 



4.1 Microlensing within the local group 

The vast majority of discrete microlensing events observed 
to date have been detected during the various monitoring 
programs (e.g. Afonso et al. 1999; Alcock et al. 2000; see 
Paczyhski 1996 for a review) directed at either the Galactic 



centre or the Magellanic Clouds, for which d os ^ 60 kpc. 
Pixel lensing studies of M 31 (e.g. Gyuk & Crotts 2000; 
Cheongho & Gould 1996) have extended the distance scale 
to ~ 1 Mpc, but this is still too nearby to be very sensitive 
to MONDian effects. Equation @ and Fig. | imply that 
9e/8m & 0.005 for all microlensing events within the local 
group. Basic lensing theory (e.g. Schneider et al. 1992) is 
then sufficient to show that only sources magnified by less 
than one per cent would be subject to any additional MON- 
Dian boost. To detect these changes, relative photometry 
accurate to ~ 0.001 mag would be required, and so there 
is little chance of constraining alternative gravity theories 
from the data currently available on single light-curves. 

Future observations of local group microlensing events 
might include very accurate photometric monitoring for long 
periods after the magnification peak, but the distinction be- 
tween MONDian and relativistic curves would only become 
apparent on scales where other deflectors along the line-of- 
sight start to have an effect. One promising avenue of inves- 
tigation is to 'stack' the light-curves of simple microlensing 
events (i.e. in which both the deflector and source are point- 
like); this is discussed in a more general context in Mortlock 
& Turner (2001). Such an approach would not only allow 
a distinction to be made between GR and MOND, but it 
would also facilitate a measurement of fyi(x) (defined in 
Section 13). 



4.2 Cosmological microlensing 

For sources at redshifts of <; 1 the Einstein radius and 
MOND radius of a point-mass deflector become comparable 
[equation (^6|) and Fig. ^| . The resultant microlensing light- 
curves are then quite different, as illustrated in Fig. ^j, and 
there is the possibility of a clean and simple test to distin- 
guish between GR, MOND, and other gravitational theories 
(Mortlock & Turner 2001). 

In order to test these predictions, accurate photome- 
try of cosmological microlensing by an isolated deflector is 
required. Such data could be forthcoming from the quasar 
monitoring programs being undertaken by Walker (1999) 
and Tadros, Warren & Hewett (2001). The idea of both 
experiments is to observe quasars seen through galactic or 
cluster haloes in the hope of a positive detection of compact 
dark matter. If MOND is correct then many of these targets 
should show no signs of microlensing, as they lie along empty 
(i.e. dark) lines-of-sight, which should be devoid of any po- 
tential deflectors in this model (e.g. Walker 1994; Mortlock 
& Turner 2001). It is thus quite likely that serendipitous 
observations of microlensing may provide the desired data, 
and indeed the spurious peak in the light-curve of gamma 
ray burst (GRB) 000301C (Sagar et al. 2000, and references 
therein) may be the result of microlensing (Garnavich, Loeb 
& Stanek 2000). Unfortunately the uncertainties in the pho- 
tometry and lack of knowledge of the source geometry pre- 
vent any discrimination between MOND and GR (Mortlock 
& Turner 2001). Nonetheless, more events should be forth- 
coming, and cosmological microlensing has the potential to 
be a very clean method by which to measure the deflection 
law of a point- mass. 
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5 QUASAR LENSING 

Since the identification of the first gravitationally-lensed ob- 
ject - Q 0957+561 (Walsh, Carswell & Weymann 1979) - 
there have been more than fifty quasar lenses discovered. 
The image configurations have been used to map the dark 
matter in the lens galaxies (e.g. Chen, Kochanek & He- 
witt 1995; Chae, Turnshek & Khersonsky 1998; Keeton, 
Kochanek & Falco 1998), and the fraction of quasars which 
are lensed places strong limits on the cosmological model 
(e.g. Kochanek 1996). In the context of a MONDian uni- 
verse, individual lenses can be used to constrain t he d eflec- 
tion law of an extended mass distribution (Section 5.1) and 
the lens stat istics are similarly revealing in a global context 
(Section 
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Figure 6. Microlcnsing light-curves, Ara(t) = mo — m(t), calcu- 
lated assuming MOND (solid lines) and GR (dashed lines). Each 
panel shows a different peak magnification (as labelled), and in 
each case the Einstein radius (i.e. the mass of the deflector) and 
minimum impact parameter are chosen to best fit the Newtonian 
result. For each value of Am ma x light-curves are for 0e/$m = 0.0 
(i.e. GR), e B /9 M = 0.1 and E /6» M = 1-0. 



5.1 Lens modeling 

If there is no dark matter then it must be possible to model 
any quasar lens with a mass distribution that, when inte- 
grated along the line-of-sight, also matches the observed sur- 
face brightness of the deflector. Other mass concentrations 
along the line-of-sight may also have some influence (e.g. 
Keeton, Kochanek & Seljak 1997), but this is potentially a 
very powerful test of MOND. 

The fact that the three-dimensional mass distribution 
is required to perform such a test suggests elliptical galax- 
ies as deflectors, and several cases of lensing by (appar- 
ently) isolated ellipticals have already been discovered: e.g. 
MG 1654+1346 (Langston et al. 1989); MG 1549+3047 
(Lehar et al. 1993); HST 1253-2914 (Ratnatunga et al. 
1995); and HST 14176+5226 (Ratnatunga et al. 1995). As- 
suming that the deflector is a triaxial ellipsoid, the con- 
straints implied by the observed surface brightness distri- 
bution leave only one free parameter, which is essentially 
the unknown extension of the galaxy along the line-of-sight. 
This is sufficiently restrictive that a successful fit would 
probably represent positive evidence in favour of MOND. 
However preliminary investigations suggest that such mod- 
eling cannot always produce a good fit to both photometric 
and lensing data, which would pose serious problems for the 
MONDian lensing theory put forward in Section ^. 

5.2 Lens statistics 

The fraction of high-redshift sources that are multiply- 
imaged increases rapidly with the cosmological constant as 
the co-moving volume element in an Qa -dominated uni- 
verse is so high (Turner 1990). Kochanek (1996) used a selec- 
tion of optical data and sophisticated modeling to show that 
the low observed number of lenses implies £Ia < 0.65 (with 
95 per cent confidence), a limit which is only marginally con- 
sistent with the conventional dark matter-dominated uni- 
verses preferred by CMB data and high-redshift supernova 
observations (e.g. Efstathiou et al. 1999). As a MONDian 
universe is likely to obey the Friedmann equations (Sanders 
1998), these lensing results are completely at odds with the 
S7 mo ~ 0.01, Qa — 0.99 model described in Section [j]. With 
the matter density limited by nucleosynthesis, the universe 
must be dominated by the vacuum energy density in order 
to satisfy the CMB power spectrum (McGaugh 2000). Thus 
it seems that no single cosmological model can explain lens 
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statistics, supernova data and CMB anisotropies within the 
framework of MOND. 

One possible flaw in the above arguments is the assumed 
form of the angular diameter distance, which is not known 
with any certainty in MOND. The lens statistics are mildly 
sensitive to the distance formula used (Ehlers & Schneider 
1986); the main aim of the supernova observations is to mea- 
sure the luminosity distance; and the angular scale of the 
first peak in the CMB power spectrum only implies a flat 
universe if the angular diameter distance (on degree scales) 
obeys the standard formulae (e.g. Carroll et al. 1992). De- 
spite this formal ambiguity, the angular diameter distance 
would have to differ quite radically from their form in GR 
to prevent the low number of quasar lenses observed being 
a serious problem for MOND. 



6 CLUSTER LENSING 

The central regions of galaxy clusters are very effective grav- 
itational lenses, as evidenced by giant arcs, the strongly 
distorted images of background galaxies (e.g. Kneib et al. 
1996). Combined with the more numerous weakly sheared 
images of background sources, the projected mass distri- 
butions of lensing clusters have been convincingly recon- 
structed in a number of cases (e.g. Mellier et al. 1993; Ham- 
mer et al. 1997; AbdelSalam et al. 1998). The inferred mass 
maps are generally correlated with the distribution of galax- 
ies within the clusters, but there is also strong evidence for 
dark matter. 

In MOND, on the other hand, the lensing signal must 
be solely due to the cluster members (and any gas present) . 
Just as lensing by isola ted galaxies offers a powerful test 
of MOND (Section 5.1), so does cluster lensing. The cor- 



relation between the mass maps and the cluster galaxies 
is certainly consistent with expectations, but the presence 
of a large number of deflectors makes a more quantitative 
analysis problematic. Most fundamentally, the net effect of 
multiple lenses is not known with any certainty, despite be- 
ing well-defined in the particular formalism described in Sec- 
tion |i| Further, even if the lens theory was known, the lack of 
a thin-lens approximation results in a mass of degeneracies, 
with the line-of-sight location of each cluster galaxy being a 
free parameter. In some cases the discrepancy between the- 
ory and observation might be so severe that this ambiguity 
is unimportant (e.g. Beckenstein & Sanders 1994), but for 
the moment no strong conclusions regarding MOND can be 
drawn from cluster lensing. 



This formalism was extended in Section g, and can be ap- 
plied to isolated mass distributions with some confidence, 
but, due to the counter-intuitive nature of MOND, cannot 
be extrapolated to multiple deflectors with any certainty. A 
further complication is the failure of the thin-lens approx- 
imation, which means that any ambiguities in the (three- 
dimensional) luminosity density of a deflector flow through 
to its lensing properties. These difficulties notwithstanding, 
this tentative theory of MONDian gravitational lensing is 
subject to a number of observational tests. 

MOND is consistent with observations of galaxy-galaxy 
lensing, although tests for deviation from azimuthal symme- 
try in the shear signal should be able to discriminate unam- 
biguously between dark matter and most alternative gravity 
theories (Section^; Mortlock & Turner 2001). Low-optical 
depth microlensing within the local group is unlikely to be a 
particularly useful in this context, but on cosmological scales 
will be a very clean probe of MOND, if sufficiently many 
lensing events are detected (Section |ij) . If MOND is consis- 
tent with such simple lensing observations, more complex 
scenarios, such as strong lensing by galaxies and clusters, 
should provide additional constraints (Sections |B| and |(| re- 
spectively). Finally, the frequency of strong lensing events is 



already at odds with MOND, as the fi n 



0.01, Q Ao ^ 0.99 



MONDian cosmology implied by CMB observations (Mc- 
Gaugh 2000) should result in many mo re multiply-imaged 
sources than are observed (Section 5.2), although there is 
some ambiguity in the redshift-distance relationship in such 
a universe. 

On balance the formalism described in Section g is a 
reasonable hypothesis for gravitational lensing within the 
framework for MOND, and must be qualitatively correct for 
isolated deflectors. The theory is subject to a number of 
tests, although most await the completion of current sur- 
veys or further theoretical development. However it is likely 
that any fully relativistic extension of MOND must be com- 
pletely non-linear to explain all of the above manifestations 
of gravitational lensing. 
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7 CONCLUSIONS 

MOND is a modification of inertia (or gravity) which can 
explain a variety of astronomical observations that would 
otherwise imply the existence of large amounts of dark mat- 
ter. Lacking a relativistic analogue, however, MOND makes 
no definite predictions about cosmology or photon propa- 
gation. Sanders (1998) took an empirical approach to the 
former problem with some success, and, in a similar spirit, 
Mortlock & Turner (2001) showed that the existing galaxy- 
galaxy lensing data are consistent with the simple MONDian 
theory of gravitational lensing described by Qin et al. (1995). 
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